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Abstract 

We define U(n) gauge theory on fuzzy 5)^ x Sjj as a multi-matrix model, which 
reduces to ordinary Yang-Mills theory on S 2 x S 2 in the commutative limit N —> oo. 
The model can be used as a regularization of gauge theory on noncommutative Rj 
in a particular scaling limit, which is studied in detail. We also find topologically 
non-trivial 17(1) solutions, which reduce to the known “fluxon” solutions in the limit 
of Rg, reproducing their full moduli space. Other solutions which can be interpreted 
as 2-dimensional branes are also found. The quantization of the model is defined 
non-perturbatively in terms of a path integral which is finite. A gauge-fixed BRST- 
invariant action is given as well. Fermions in the fundamental representation of the 
gauge group are included using a formulation based on 50(6), by defining a fuzzy 
Dirac operator which reduces to the standard Dirac operator on 5 2 x 5 2 in the 
commutative limit. The chirality operator and Weyl spinors are also introduced. 
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1 Introduction 


Gauge theories on noncommutative spaces have received much attention in recent years. 
One of the reasons is the natural realization of such theories in the framework of string 
theory and D-branes [1], however they deserve interest also in their own right; see [2,3] 
for some reviews. One of the most remarkable new features of noncommutative gauge 
theories is the fact that they can be defined in terms of multi-matrix models, which means 
that the action involves only products of “covariant coordinates” X, = Xi + Ai, with gauge 
transformations acting as X{ —> UXiU -1 . In particular for certain quantized compact 
spaces such as fuzzy spheres and tori, these Xi are finite-dimensional Hermitian matrices 
of size N. Nevertheless, the conventional gauge theory is correctly reproduced in the limit 
N —> oo. This leads to a natural quantization prescription by simply integrating over 
these matrices. For the much-studied case of the quantum plane WL the matrices Xi are 
infinite-dimensional, and the precise definition of the models is quite non-trivial. This is 
particularly obvious by noting that the naive action for gauge theory on Mg contains sectors 
with any rank of the gauge group U (n) [4]. To have a well-defined theory and quantization 
prescription, a regularization of gauge theory on Mg based on the finite compact case is 
therefore very desirable. Furthermore, the formulation as multi-matrix model leads to 
the hope that non-trivial results may be obtained using the sophisticated techniques from 
random matrix theory. We introduce in this paper such a matrix model for fuzzy S 2 X S 2 , 
and study its relationship with Mg. 

In the 2-dimensional case, this matrix-model approach to gauge theory has been studied 
in considerable detail for the fuzzy sphere [5-10] and the noncommutative torus Tg 
[11-14], both on the classical and quantized level. It is well-known that Mg can be obtained 
as scaling limits of these spaces S]^ and T 2 r at least locally, which suggests a correspondence 
also for the gauge theories. This correspondence of gauge theories has been studied in great 
detail for the case of Hg —> Mg [12,15,16] on the quantized level, exhibiting the role of 
certain instanton contributions. A matching of gauge theory on the classical level can also 
be seen for —► Mg [17, 18], which is implicitly contained in Section Q of the present 

paper. 

In 4 dimensions, the quantization of gauge theory is more difficult, and a regularization 
using finite-dimensional matrix models is particularly important. The most obvious 4- 
dimensional spaces suitable for this purpose are T , S' 2 x S 2 and CP 2 . On fuzzy CP 2 r 
[19-21], such a formulation of gauge theory was given in [22], This can indeed be used to 
obtain Mg for the case of U( 2) -invariant 9 VJ . The case of M 2 x S ]y as regularization of Mg 
with degenerate Qjj was considered in [18,23], exhibiting a relation with a conventional 
non-linear sigma model. A formulation of lattice gauge theory for even-dimensional tori 
has been discussed in [14,24,25]. Related “fuzzy” solutions of the string-theoretical matrix 
models [26] were studied e.g. in [27,28], see also [29]. 

In the present paper we give a definition of U (n) gauge theory on fuzzy x S'fr , which 
can be used to obtain any Mg as a scaling limit. The action is a simple generalization of the 
matrix model approach of [7] for fuzzy . It differs from similar string-theoretical matrix 
models [26] by adding a constraint-term, which ensures that the “vacuum” solution is 
stable and describes the product of 2 spheres. The fluctuations of the covariant coordinates 
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then correspond as usual to the gauge fields, and the action reduces to ordinary Yang-Mills 
theory on 5 2 x 5 2 in the limit N —* oo. The quantization of the model is defined by a 
finite integral over the matrix degrees of freedom, which is shown to be convergent due to 
the constraint term. We also give a gauge-fixed action with BRST symmetry. 

We then discuss some features of the model, in particular a hidden 50(6) invariance 
of the action which is broken explicitly by the constraint. This suggests some alternative 
formulations in terms of “collective matrices”, which are assembled from the individual 
covariant coordinates (matrices). This turns out to be very useful to construct a Dirac 
operator, and may help to eventually study the quantization of the model. The stability of 
the model without constraint is also discussed, and we show that the only flat directions of 
the 50(6) -invariant action are fluctuations of the constant radial modes of the 2 spheres. 

As a further test of the proposed gauge theory, we study in Section |H1 topologically 
non-trivial solutions (instantons) on 5]y x Sfj. We find in particular a simple class of 
solutions which can be interpreted as 0(1) instantons with quantized flux, combined with 
a singular, localized “flux tube”. They are related to the so-called “fluxon” solutions of 
0(1) gauge theory on Mg. Solutions which can be interpreted as 2-dinrensional spherical 
branes wrapping one of the two spheres are also found and are matched with similar 
solutions on Mg. We then study the relation of the model on 5^ x 5^ with Yang-Mills 
theory on Mg, and demonstrate that the usual Yang-Mills action on Mg is recovered in the 
appropriate scaling limit. Some aspects of 0(1) instantons (“fluxons”) on Mg are recalled 
in Section o and we show in detail how they arise as limits of the above non-trivial 
solutions on 5^ x 5^-. In particular, we are able to match the moduli space of n fluxons, 
corresponding to their location on Mg resp. 5)y x 5)y. We find in particular that even 
though the field strength in the “bulk” vanishes in the limit of Mg, it does contribute 
to the action on 5^ x Sfj with equal weight as the localized flux tube. This can be 
interpreted on Mg as a topological or surface term at infinity. Another unexpected feature 
on Sjj x Sff is the appearance of certain “superselection rules”, restricting the possible 
instanton numbers. In other words, not all instanton numbers on Mg are reproduced for a 
given matrix size JV, however they can be found by considering matrices of different size. 
This depends on the precise form of the constraint term in the action, which is hence seen 
to imply also certain topological constraints. To recover the full space of ADHM solutions 
on Mg starting from 5^ x 5^ remains an open challenge, which is non-trivial since the 
concept of self-duality does not extend naturally to the fuzzy case. 

We should mention here that topologically non-trivial configurations have also been 
discussed more abstractly in terms of projective modules using a somewhat different for¬ 
mulation of gauge theory on fuzzy spaces, see in particular [30,31]. 

Finally in Section [H] we include charged fermions in the fundamental representation of 
the gauge group, by giving a Dirac operator D which in the large N limit reduces to the 
ordinary gauged Dirac operator on 5 2 x 5 2 . This Dirac operator is covariant under the 
50(6) symmetry of the embedding space 5 2 x 5 2 C K 6 , and exactly anti-commutes with 
a chirality operator. The 4-dimensional physical Dirac spinors are obtained by suitable 
projections from 8-dimensional 50(6) spinors. This projection however commutes with 
D only in the large N limit, and is achieved by giving one of the 2 spinors a large mass. 
Weyl spinors can then be defined using the exact chirality operator. 
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2 The fuzzy spaces S% and S% L x Sjf 


We start by recalling the definition of the fuzzy sphere in order to fix our conventions and 
notation. The algebra of functions on the fuzzy sphere is the finite algebra S^ generated 
by Hermitian operators X{ = (x\,X 2 ,X 3 ) satisfying the defining relations 


\xi-, Xj\ iA-jyCijkXkj ( 1 ) 

x\ + x\ + x 2 = R 2 . (2) 


They are obtained from the iV-dimensional representation of su(2) with generators A* (i = 
1,2,3) and commutation relations 


[Aj, Xj] — Rijk^ki 


3 

A?;Aj 

2—1 



(see Appendix El by identifying 

Xi — A^v A i, 


A N 


2 R 

VN 2 - 1' 


(3) 

(4) 


The noncommutativity parameter A^v is of dimension length. The algebra of functions Sj (r 
therefore coincides with the simple matrix algebra Mat(N, C). The normalized integral 
of a function / E 5^ is given by the trace 


/ 


/ 


A-kR? 

N 


tr (/)■ 


(5) 


The functions on the fuzzy sphere can be mapped to functions on the commutative sphere 
S 2 using the decomposition into harmonics under the action 


Jif = [Kf] 


( 6 ) 


of the rotation group 517(2). One obtains analogs of the spherical harmonics up to a 
maximal angular momentum IV — 1. Therefore 5^ is a regularization of S 2 with a UV 
cutoff, and the commutative sphere 5 2 is recovered in the limit N —> oo. Note also 
that for the standard representation CUB. entries in the upper-left block of the matrices 
correspond to functions localized at X 3 = R. In particular, the fuzzy delta-function at the 
“north pole” is given by a suitably normalized projector of rank 1 , 



N JV-1 JV-1 
AnR 2 2 2 


(7) 


where \ N 0 1 ) is the highest weight state with maximal eigenvalue of A 3 . Delta-functions 
with arbitrary localization are obtained by rotating 0 . 
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The simplest 4-dimensional generalization of the above is the product S'fr^ x S^ R of 
2 such fuzzy spheres, with generally independent parameters N^.r. It is generated by a 
double set of representations of su( 2) commuting with each other, i. e. by Af, X f satisfying 


[Af.A f] = ie ijk Af, [Af, Xf ] = ie ijk X k , 

[Af,A f] = 0 

for i, j = 1, 2, 3, and Casimirs 




L_Nl-l 


i= 1 


B N 2 - 1 
J^AfAf = -L. 

1=1 


This can be realized as a tensor product of 2 fuzzy sphere algebras 


( 8 ) 


Af 

= Xi <g> 1n r xn r , 

(9) 

Af 

= 1n l xN l <8> Aj, 

(10) 

X S Nr 

= Mat (A/, C) where 



Af = N l N r . 

(ID 


The normalized coordinate functions are given by 


l,r 21? LR 

X,- = — / - A,- , 

* v/(ATL,R)2 _ 1 * 


I>.T = -R 2 = E<T) 2 - 


( 12 ) 


This space 1 can be viewed as regularization of S 2 x 5 2 C M 6 , and admits the symmetry 
group SU(2)l x SU(2) r C 50(6). The generators xf’ R should be viewed as coordinates 
in an embedding space R 6 . The normalized integral of a function / € S^ L x 5]y is now 
given by 



167r 2 i? 4 

Af 


tr (/) 


V_ 

Af 


t r(/), 


(13) 


c2 v c2 

AT ^ AT- 


where we define the volume V := 167r 2 i? 4 . We will mainly consider Nr = Nr in the 
following. 


2.1 The quantum plane limit Rg 

It is well-known [32] that if a fuzzy sphere is blown up near a given point, it can be used 
to obtain a (compactified) quantum plane: Consider the tangential coordinates x\^ near 
the “north pole”. Setting 

R 2 = Nd/2, (14) 

1 In principle one could also introduce different radii R l,r for the 2 spheres, but for simplicity we will 
keep only one scale parameter R (and usually we will set R = 1). 
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they satisfy the commutation relations 


[x\,x 2 ] = \Jr 2 - x\ - x\ = i 6 (l + 0(l/N)). (15) 

Therefore in the large N limit with m keeping 6 fixed, we recover 2 the commutation 
relation of the quantum plane, 

[xi,x 2 ] = iO (16) 

up to corrections of order jj. Similarly, starting with Sfj L x Sjf and setting 

R 2 = N L , R e LjR / 2, (17) 

we obtain in the large N R , N r limit 

[xf,xf] = ieijd L , [Xj t ,xf]=ie i jO R , 

[xi,xf\ = 0. (18) 

This is the most general form of with coordinates (xi,..., X 4 ) = (xf, x|s x R , x R ) (after 
a suitable orthogonal transformation). The integral of a function fix) then becomes 

J fix) -> 4vr 2 6' L 6 , i?tr(/(x)) =: J /(x), (19) 

Sjf L X S% R Rg 

which has indeed the standard normalization, giving each “Planck cell” the appropriate 
volume. 


3 Gauge theory on fuzzy S 2 x S' 2 


We start with the most general case, and construct a matrix model having S'j^ L x Sjj R as 
its ground state. The fluctuations around this ground state will produce a gauge theory. A 
simplified and more elegant formulation in terms of “collective matrices” similar as in [7] 
for the fuzzy sphere will be given later in Section 0J 

In the fuzzy case, it is natural to construct S' 2 x 5^ as “submanifold” of M 6 . We 
therefore consider a multi-matrix model with 6 dynamical fields (“covariant coordinates”) 
Bj t J and B R (i = 1,2,3), which are J\f x J\f Hermitian matrices. As action we choose the 
following generalization of the action in [7,8], 


5 = 


cy Ra jbRia jb + 7T + VTfi 


( 20 ) 


with a,b = L, R and i.j = 1,2,3; summation over repeated indices is implied. Here (Pl,r 
are defined as 


._ 1 (pLpL 

n - 7 B, - 


.. ._ 1 (pRpR N R~ 1 \ 

VR — ^2 i B i B i - 4 - h 


( 21 ) 


2 One could be more sophisticated and use the stereograpliic projections as in [32], which leads essentially 

to the same results. 
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and the terms </? 2 + tp\ in the action ensure that the unwanted radial degrees of freedom 
are suppressed [7,8]. R denotes the radius of the two spheres, which we keep explicitly to 
have the correct dimensions. The held strength is defined by 

FiLjL = ] + EijkB k ), 

FiRjR = ,Bj ] + tijkBfc ), 

FiLjR = ^{i[B R ,B R }). (22) 

This model m is manifestly invariant under SU[2)l x SU(2)r rotations acting in the 
obvious way, and U(M) gauge transformations acting as B R,R — * UB R,R U~ 1 . We will see 
below that this reduces indeed to the U (1) Yang-Mills action on S 2 x S 2 in the commutative 
limit. Note that if the action (120)1 is considered as a matrix model, the radius drops out 
using m- The equations of motion (e.o.m.) for B^ are 

/y 2 — 1 

{Bt B R Bf - } + (B[ + it iji,. I>j />’/:) 

+Rijk[B R , (B% + ie^rsB R B R )] + [B R , [B R , B R ]\ = 0, (23) 

and those for B^ are obtained by exchanging L R. By construction, the minimum or 
ground state of the action is given by F = tp = 0, hence B i ' = A- ’ as in (J9J, (flOl) up to 
gauge transformations; cp. [22] for a similar approach on CP 2 . We can therefore expand 
the “covariant coordinates” Bf and Bf- around the ground state 


Bf = A f + RAf, 


(24) 


A R ’ R -4 A ,Rli = UA^u- 1 + U[Xf’ H , U~\ 


^/L,R 
i 


L R 

where a G {L, R} and Af is small. Then A- ’ transforms under gauge transformations as 

(25) 

and the held strength takes a more familiar form 3 , 

F,LjL = + 


R 

FiRjR = i([f.Yl-[f.Y] + [-4f,Yl). 

FiRiR = i([f,Yl-[f.Y] + !YY])- 


(26) 


So far, the spheres are described in terms of 3 Cartesian covariant coordinates each. In the 
commutative limit, we can separate the radial and tangential degrees of freedom. There 
are many ways to do this; perhaps the most elegant for the present purpose is to note 

3 We do not distinguish between upper and lower indices L, R. 
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that the terms J cp 2 L + ip 2 i in the action imply that pl,r is bounded for configurations with 
finite action. Using 


\ L 

^ 


\L 


L aL 


+ Af-^ + A?A 


(27) 


and similarly for pp it follows that 


XiA? + AiXi 


o( 


N> 


(28) 


for finite A?. This means that Af is tangential in the (commutative) large N limit. 
Alternatively, one could consider (f>L = Nipp, which would acquire a mass of order N 
and decouple from the other fields 4 . The commutative limit of (1201) therefore gives the 
standard action for electrodynamics on S 2 x S' 2 , 

c _ _ f ppt 

° 2 g 2 J iajb r iajb 

S 2 xS 2 


with a,b = L,R and i,j = 1, 2,3. Here F? L - R denotes the usual tangential field strength. 

L R 

This can be seen most easily noting that e.g. at the north pole x 3 ’ = R, one can replace 


,L,R 

\ _ 

R ' 


d 


—£■ 


' 3 th*' R 


(29) 


in the commutative limit, so that upon identifying the commutative gauge fields A, 

A (cl)L,R _ .L,R 


(cl) 


via 

(30) 


the field strength is given by the standard expression F- L - R = d L A {d)R _ QR A (cl)L 

etc. 


U(k) gauge theory 

The above action generalizes immediately to the nonabelian case, keeping precisely the 
same action m , <ED but replacing the matrices B^’ R by kAf x kAf matrices, cp. [7]. 
Expanding them in terms of (generalized) Gell-Mann matrices, the same action (1201) is the 
fuzzy version of nonabelian U(k) Yang-Mills on S 2 x S' 2 . 


4 A formulation based on 50(6) 

The above action can be cast into a nicer form by assembling the matrices B^’ R into 
bigger “collective matrices”, following [7]. Since it is natural from the fuzzy point of 
view to embed S 2 x S 2 C i? 6 with corresponding embedding of the symmetry group 
SO(3)l x SO(3)r C 50(6), we consider 

_ B, = (B\\ B r ) (31) 

4 The constraints ifiL = 0 = <pu could also be imposed by hand; however the suppression through the 
above terms in the action is more flexible, as we will see in Section 0 
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(Greek indices //, v denoting from now on all the six dimensions) to be the 6 -dimensional 
irrep of so( 6 ) = su( 4). Since (4)<8>(4) = (6)©(10), it is natural to introduce the intertwiners 

7m = (7* L ,7* R ) = ( 32 ) 

of ( 6 ) C (4)<8> (4). where a, f3 denote indices of (4). We could then assemble our dynamical 
fields into a single AJ\[ x 4A7 matrix 


B = Bfj 7 ^ + const • 1. (33) 

Of course the most general such 4A f x 4A7 matrix contains far too many degrees of freedom, 
and we have to constrain these B further. Since SU(4) acts on B as B —> U T BU, 
the 7 ^ can be chosen as totally anti-symmetric matrices, which precisely singles out the 


(6) C (4) <S) (4). One can moreover impose 


^7 

t-i 

II 

II 

1 

(34) 

and 


7 i 7 j = $ij T i^ijkIk 1 

(35) 

1 *1? = -Sij-Cijklk, 

(36) 

O 

II 

(37) 


which will be assumed from now on; we will give two explicit such representations in cm 
em This would suggest to constrain B to be antisymmetric. However, the component 
fields Bp are naturally considered as Hermitian rather than symmetric matrices. Fur¬ 
thermore, since the 7 ^ = ( 7 ^)“have two upper indices, they do not form an algebra. 
There are now 2 ways to proceed. We can either separate them again by introducing two 
4A f x 4A7 matrices, 

B L = \ + Bht, B R = - + B R 1 R (38) 

breaking 50(6) —► 50(3) x 50(3). This will be pursued in Appendix [BJ Alternatively, 
we can use the 7 ^ with the above properties to construct the 8 x 8 Gamma-matrices 


pM 


0 7^ \ 

7^t 0 ) ’ 


which generate the 50(6)-Clifford algebra 


{T M , T"} 


y/My/O _j_ y/^y/Mt Q 

Q ryiA ry V _|_ ry V t/-yM 


2 8 ^. 


This suggests to consider the single Hermitian 8 A f x 8 A f matrix 


o = + Oo 


0 b l \ ( 0 b r 

B l 0 ) + { —B r 0 


= : C L + C R , 


(39) 


(40) 
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(41) 


where Co = Cq + C R denote the constant 8 x 8-matrices 

= 4 rfr ^ = Kl o)’ (42) 

c 0 R = -^rfr«rf = ^_° 1 J) (43) 

in the above basis. This is very close to the approach in [7], and using the Clifford algebra 
and the above definitions one finds indeed 

C 2 = + i + ^F IW . (44) 

Here £g^ = — |[r^,r jy ], and the field strength F )W coincides with the definition in (1221) if 
written in the L — R notation, 

B{a jb — i\Biai B T Sab^-ijk B } :a . 


Therefore the action 

N 2 N 2 — 1 

5 6 = Tr((C 2 - —) 2 ) = Str^i^- —f + 4tr F^F^ (45) 

is quite close to what we want. The only difference is the term (B^B^ — Af ' 2 ~ 1 ) 2 instead 

of (B iL B iL - + ( B iR B iR - ) 2 , for 21V 2 = JV| + 1V|. This difference is easy 

to understand: since (ESI is 50(6)-invariant, the ground state should be some 5 5 . We 
therefore have to break this SO ( 6 )- invariance explicitly, which will be done in the next 
section. However before doing that, let us try to understand action ESI) better and see 
whether it leads to a meaningful 4-dimensional field theory. We show in Appendix □ 
by carefully integrating out the scalar components of B R)R that the 50(6)- invariant 
constraint term in m induces the second term in the following effective action 

Sq S ~ 4tr ( F tw F, w - ( F iL Xi L - F iR x iR )—j— 1 ( F iL x iL - F iR x iR ) ) (46) 

V 4 v2 d A) J 

in the commutative limit, where Fn = \€ijkFj R kL etc. Comparing the second term with 
we see that the zero mode of the Laplace operator can produce a contribution 
that cancels the corresponding contribution from F^ U F^ U , but that all higher modes are 
smaller by at least a factor of 2 ( 4 — d^d^). Therefore, the action (HH1) is positive definite 
except for the obvious zero mode 5B= e, SB? 1 = —e. This means that the geometry 
of S\ x S 2 r is locally stable even with the S O ( 6 )-symmetry unbroken, except for opposite 
fluctuations of the radii. 

4.1 Breaking 50(6) —> 50(3) x 50(3) 

To obtain the original action m f° r s 2 X S 2 , we have to break the 50(6)-symmetry down 
to 50(3) x 50(3). We can do this by using the left and right gauge fields C L and C R 
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introduced in m separately. Their squares are 




Bn Bn + — + 

BiRBiR + - - 



As both 7 l L , 7 ^ and Yl'Yr are traceless, we have 



With these terms we can recover our action as 



which is precisely the action CT for gauge theory on S]^ L x omitting the overall con¬ 
stants. Hence the action is formulated as 2-matrix model, however with highly constrained 
matrices Cl,Cr. This formulation using the Gamma-matrices is very natural and useful 
if one wants to couple the gauge fields to fermions, as discussed in Section [HJ 
For simplicity, we will only consider Nr = Nr = N from now on. 


5 Quantization 

The quantization of the gauge theory defined by m or its reformulation (EH is straight¬ 
forward in principle, by a “path integral” over the Hermitian matrices 



(48) 


Note that there is no need to fix the gauge since the gauge group U(J\f) is compact. The 
above path integral is well-defined and finite for any fixed J\f. To see this, it is enough 
to show that the integral J dB^ exp(—( B- J — (N^ — l)/4 ) 2 — (B^B^ — ( N 2 — l)/4) 2 ) 
converges, since the contributions from the field strength further suppress the integrand. 
This integral is obviously convergent for any fixed N. 

For perturbative computations it is necessary to fix the gauge, and to substitute gauge 
invariance by BRST-invariance. Such a gauge-fixed action will be presented next. 

5.1 BRST Symmetry 

To construct a gauge-fixed BRST-invariant action, we have to introduce ghost fields c 
and anti-ghost fields c. These are fermionic fields, more precisely A f x N— matrices with 
entries which are Grassman variables. 
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The full gauge-fixed action reads: 


Sbrst = S+ -^=tr(c[A M , [Bp,c]] - (^b - [A M , B^\)b ), 

where b is an auxiliary (Nakanishi-Lautrup) field. This action is invariant with respect to 
the following BRST-transformations: 

sB^ = [Bf, , c] sc = cc (49) 

sc = b sb = 0 (50) 

(matrix product is understood), where the BRST-differential s acts on a product of fields 
as follows: 

s(XY) = X(sY) + {-l) eY (sX)Y. 

Here ey denotes the Grassman-parity of Y 

_ ( 0 Y bosonic 
65 [ 1 7 fermionic. 

As usual, it is not difficult to check that these BRST-transformations are indeed nilpotent, 
i.e. 

s 2 = 0. 

Integrating out the auxiliary field b leads to the following action 

‘S'brst = S + ■^tr(c[A #1 , [Bp, c]] — — [A M , B^\[ A„, B u ]). 

Setting a = 1 corresponds to the Feynman gauge. This is indeed what one would obtain 
by the Faddeev-Popov procedure. The action S' is invariant with respect to the following 
operations: 


s'B^ = [Bp, c] 
s'c = cc 
sc — [ , Bp\. 

Since we have used the equations of motion of b, the BRST-differential s' is not nilpotent 
off-shell anymore but still we have 

„/2 I _ rj 

* Ion-shell -u ' 

6 Topologically non-trivial solutions on Sjj x Sjj 

In order to understand better the non-trivial solutions found below, we first note that the 
classical space S 2 x S 2 is symplectic with symplectic form 

u> = uj l + lo r , ( 51 ) 
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where 


( 52 ) 


LO j ^ jyi^ ( Jjk- 1 ' 1 dXj dx k 
and similarly lo r . The normalization is chosen such that 

[ U L ’ R = 1= [ u L A u> R ( 53 ) 

JS\ R Js 2 xS 2 

so that lu l ,lu r generate the integer cohomology H*(S 2 x S 2 ,Z). Noting that lo is self¬ 
dual while Co := oj l — uj r is anti-selfdual, it follows immediately that both F = 2-koj and 
F = 27 too are solutions of the Abelian field equations. More generally, any 

p(m L ,m R ) _ 2TrniLUJ L + 2nmRLO R (54) 

for any integers mt, vir is a solution. In bundle language, they correspond to products of 
2 monopole bundles with connections and monopole number rriL,R over S\ R . Following 
the literature we will denote any such non-trivial solution as instanton. 


6.1 Instantons and fluxons 


We are interested in similar non-trivial solutions of the e.o.m. (El in the fuzzy case. The 
monopole solutions on the fuzzy sphere S 2 N are given by representations of su( 2 ) 

of size N — m [33], which lead to the classical monopole gauge fields in the commutative 
limit as shown in [7]. It is hence easy to guess that we will obtain solutions on Sfj x Sjf 
by taking products of these: 

B r = a L \*~ mL ( 55 ) 

B r = a R t N . mL ®\^- mR ( 56 ) 

where A ^ mL ’ R are the N — m,L,R dimensional generators of su( 2). It is not difficult to 
verify that these are solutions of (El with a L,R = 1 + for rriL,R N, with field 
strength 

m L L m R r 

FiLjL — ~ 2i?3 5 FiRjR = ^ijk^k i ^iLjR = 0? (57) 

while B • B — jV 4 ~ 1 —> 0 as N —> oo. This means that F = —2irm L uj L — 2irm R LU R in the 

commutative limit, so that indeed 


[ L = _ m FR' 

J 2 7T 

qL ,r 


( 58 ) 


Notice that the Ansatz (El implies that all matrices have size A f = (N — 171l)(N — mR ), 
which is inconsistent if we require that A f = N 2 in order to have the original Sfj x 
vacuum. Therefore it appears that these solutions live in a different configuration space, 
similar as the commutative monopoles which live on different bundles. However, the 
situation is in fact more interesting: the above solutions can be embedded in the same 
configuration spaces of N 2 x N 2 matrices as the vacuum solution if we combine them with 
other solutions, which have finite action in four dimensions 5 . They are in fact crucial to 

5 as opposed to 2 dimensions, which is the reason why they were not considered in [7] 
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recover some of the known U( 1) instantons in the limit S jy —> Mg resp. Sjf x Sjy —> M|, 
as we will see. Consider the following Ansatz 


B L ’ R 


diag(d 


L,R 
i, 1 j 



in terms of diagonal matrices (ignoring the size of the matrices for the moment). These 
are solutions of CSl in two cases, 


( j\j2 _ o . 

typeA 


0 , 


type B 


(i.e. = 0 in type B). The associated field strength is 

FiLjL = diag(4,i, dfc,n)> F LR = 0 , 


(59) 


(60) 


and a similar formula for FmjR. The constraint term is then (B ■ B — N 4 * ) —> — i for 
type A, and (B ■ B — jV 4 ~ 1 ) —> — jV 4 ~ 1 for type B in the large IV limit. In particular, only 
the type A solutions will have a finite contribution 




uxon 


V f n 2n N 2 — 3\ 8-7T 2 

4 J ~gF 


(61) 


to the action 6 , which for IV —*• oo is only due to the field strength. We will see below that 
these type A solutions can be interpreted as a localized flux or vortex, and we will call 
them “fluxons” since they will reduce in a certain scaling limit to solutions on Mg which 
are sometimes denoted as such [34-36]. 

One can now combine these “fluxon” solutions with the monopole solutions (j56l) in the 
form 


II 

cq 

( a L ® 

0 


V o 

dia g( d O’-- 

•» d £n) 

B r = 


0 


V o 

diag(d^f,.. 

■i d i,n) 


(62) 


These are now matrices of size A f = (IV — mi)(N — mu) + n, which must agree with 
J\f = N 2 , say. This is clearly possible for 


rriL = —tur = m, n = m 2 , 


(63) 


while for mi, ^ —mu the contribution from the fluxons would be infinite since n = 0(N). 
To understand these solutions, we can compute the gauge field from ( 01 ) . 

Ai = ( B i - x i ® M = A i (® L > xR )- ( 64 ) 

6 A finite action can also be obtained for the type B solution using a slightly modified action , as 
discussed below. 
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To evaluate this, we first have to choose a gauge, i.e. a unitary transformation U for (PI) 
which allows to express e.g. \^~ mL <g>l N _ mR in terms of xf oc A,^® Hat and x R oc Hat® . 
For example, in the case mi = — rir = m this can be done using a unitary map 

U : C N ~ m ® C N+m ® C m2 ® C^, (65) 


mapping a (N — m ) x (N + m) matrix into a N x N matrix by trivially matching the 

2 

upper-left corner in the obvious way, and fitting C m into the remaining lower-right corner. 
With this being understood, one can write 


RAi (x L , x R ) = (a L Af m - Af) ® 1 N+m + Af ® (1at+ to - lkv) + (d - terms) 

= A[ m \x L ) + sing(rr 3 = —R,x r = —R) ( 66 ) 


where A^ n \x L ) is indeed the gauge field of a monopole with charge m on S\ in the large 
N limit, as was checked explicitly in [7]. Here sing(x 3 = —R,x r = —R) indicates a 
field localized at the “south pole” of S\ and/or S ^ which becomes singular for large N. 
It originates both from “cutting and pasting” the bottom and right border of the above 
matrices using U (leading to singular gauge fields but regular field strength at the south 
poles), as well as the d-block (leading to a singular field strength). To see this recall that 
in general for the standard representation E3J of fuzzy spheres, entries in the lower-right 
block of the matrices correspond to functions localized at x '3 = —R, cp. 0 . The gauge 
field near this singularity will be studied in more detail in Section f7..'->l The field strength 
is 


FiLjL = ~ 


rri 


2R3 €i J kX k + e ijk £,2 


E<* 

i =1 


'L p 

k,i r i 


(67) 


in the commutative limit, where Pi are projectors in the algebra of functions on S ^ x of 
rank 1 ; recalling 0 , they should be interpreted as delta-functions Pi = jj? d^ 4 ^(x 3 = — R). 
Similar formulae hold for Af(x L ,x R ) and FiRjR, while Flr = 0. 

We assumed above that these delta-functions are localized at the south poles x% = 
x R = —R. However, the location of these delta-functions can be chosen freely using gauge 
transformations. This can be seen by applying suitable successive gauge transformations 
using N — fc-dimensional irreps of SU( 2) for k = 0,1, ...,m — 1, which from the classical 
point of view all correspond to global rotations, successively moving the individual delta- 
peaks. Therefore the solution m should in general be interpreted as monopole on S 2 x S 2 
with monopole number mi = —mR = m , combined with a localized singular field strength 
characterized by its position and a vector d^. We will see in Section 0 that it becomes 
the “fluxon” solution in the planar limit we therefore also call it a “fluxon”. 

The total action of these solutions m is the sum of the contributions from the 
monopole field plus the contribution from the fluxons m , which both give the same 
contribution 

47T^ 

S {m) = (2m 2 + 2m 2 ) (68) 

in the large N limit, using (PI). The first term is due to the global monopole field 10 ), 
and the second term is the contribution of the fluxons through the localized field strength. 
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The interpretation of these solutions depends on the scaling limit N —> oo which we 
want to consider. We have seen that in the commutative limit keeping R = const, these 
solutions become commutative monopoles on S 2 x S 2 with magnetic charges uir = —thr, 
plus additional localized “fluxon” degrees of freedom. For large R, the field strength of 
the monopoles vanishes, leaving only the localized fluxons. In particular, we will see in 
the following section that in the scaling limit x Sj^ —► Mg only the fluxons survive and 
become well-known solutions for gauge theory on Mg. Away from this localized fluxon the 
gauge field becomes a flat connection, which is however topologically nontrivial. This is 
very interesting as it shows that one can indeed use these fuzzy spaces as regularization 
for gauge theory on Mg n . 

A final remark is in order: if we fix the size J\f of the matrices, only certain fluxon and 
monopole numbers are allowed, given by Otherwise the number n of fluxons and 

hence the action would diverge with N. This can be seen as an interesting feature of our 
model: viewed as a regularization of gauge theory on Mg, this points to possible subtleties 
of defining the admissible field configurations in infinite-dimensional Hilbert spaces and 
relations with topological terms in the action. On the other hand, we could accommodate 
the most general solutions including also type B solutions m by modifying the action 
similar as in [7]. For example, 


S = 


ABpBp , T T N? — l.o 

N 2Ri ~ 4 > 


(p R p R Nr 1 2 

H-^-; 


i ? 4 


) T c ^F'ia,jbF'ia,jb 


(69) 


leads to the same commutative action, but with a vanishing action for the Dirac string in 
the type B solutions. 


6.2 Spherical branes 

Consider the following solutions 


Bp = 


L \N—m 


a L A; 


0 


0 diag 

BP = 1 N ®\? 


® lbv, 


which are matrices of size M = N 2 . The corresponding field strength is 


F, 


m 


iLjL 


2 R 3 

Frr = F lr = 0 


^ijk^k jy) P \ dk^Pi 


i =1 


(70) 


(71) 


where Pi are projectors in the algebra of functions on S\ of rank 1 which should be 
interpreted as delta-functions Pi = 4 ^ 2 = —R). In particular the gauge field A 

vanishes on Sp, while on S\ there is a monopole field together with a singularity at a 
point. This is similar to the fluxons on the previous section, but now only on S'"). This 
leads to the interpretation as 2-dimensional brane wrapping on S^, located at a point on 
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S £. The action for these solutions is infinite. In the limit x S^y —> Mg, the flux will 
be located at a 2-dimensional hyperplane. Such solutions for gauge theory on Mg were 
found in [4,37], which would be recovered in the scaling limit Sjj x 5^ —> Mg as discussed 
in Section Q In a similar way, we can interpret solutions with any as branes 

wrapping on S\ and S‘f { . 

7 Gauge theory on from S% L x S 2 Nr 

We saw in Section l2~Tl that Mg can be obtained as a scaling limit of fuzzy S^ L x Sfj . 
Here we will extend this scaling also to the covariant coordinates B^, thereby relating the 
gauge theory on S^ L x S^ R to that on Mg and hence providing a regularization for the 
latter. We will in particular relate the instanton solutions on these two spaces. 

On noncommutative Mg, all (7(l)-instantons were constructed and classified in [4]. 
They can be interpreted as localized flux solutions, sometimes called fluxons. One can 
indeed recover these instantons from corresponding solutions on Sjj, as we will show below. 
However since we are mainly interested in the 4-dimensional case here, we will only present 
the corresponding constructions on S^ L x S^ R resp. Mg here, without discussing the 2 - 
dimensional case separately. It can be recovered in an obvious way from the considerations 
below. 

The situation on Mg is more complicated, and there are different types of non-trivial 
17(1) “instanton” solutions on Mg. Assuming that 6 ^ is self-dual, there are two types 
of instantons: first, there exist straightforward generalizations of the localized “fluxon” 
solutions with self-dual field strength. These will be discussed in detail here, and we 
will show how these solutions can be recovered as scaling limits of the solutions (El on 
S^ L x Sjj R . This is one of the main results of the present paper. In particular, the moduli 
of the fluxon solutions on Mg will be related to the free parameters d^ )R in (1(1211 . This 
supports our suggestion to use gauge theory on S^ L x S%j as a regularization for gauge 
theory on Mg. However there are other types of (7(1) instantons on Mg which were found 
through a noncommutative version of the ADHM equations [38-43], in particular anti¬ 
selfdual instantons which are much less localized than the fluxon solutions. To find the 
corresponding solutions on S^ L X Sfj is an interesting open challenge. 

7.1 The action 

The most general noncommutative Mg is generated by the coordinates subject to the 
commutation relations 

[S'/zj^'i/] — i (72) 

where /i, v £ {1,..., 4}. Using suitable rotations, 9^ can always be cast in the following 
form: 

0 0i2 0 0 \ 

-012 0 0 0 

0 0 0 0 34 

0 0 -034 0 
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We will assume that @12 > 0 and $34 > 0 for simplicity in this section. Then define 

Xl ' 2 := (73) 

/20 34 (7A -. 

(75) 


^3,4 := 


^ := 83“ - ^ + T^((Bf S ) 2 + (Bf R ) 2 ) . 


which should be interpreted as a blow-up near the north pole. In the scaling limit (El, 


1 


1 


R z = -xNlOu = -N r 9 12 


00 


(76) 


the X will become the covariant coordinates on the “tangential” Mg as Nl r —> 00 , and 0 
remains an auxiliary field. To see this, we compute for the field strength 

±([B^B?]) = -l—[ Xl ,X z \, etc., 
ft (712^34 




Tinv {Xl + ^^ - ojr (Xl)2 ] } 

vi2t734Bl ^12 

(v 2 + »[w,^]-^-[W,(v 2 ) 2 ]) 

i7l2t734Bl Zt /12 


ft ( 742(734 


1 - —■ - ^ 12^34 j ^12^34 (( v 42 


(0 12 +i[*4,* 2 ] + 


-</>L - 


i ? 2 Y 2 i ? 4 

Analogous expressions hold for . For the potential term we get 


((A4) 2 + (X 2 ) 2 )). 


1 (pLpL N L 1 

J?<B- B ‘ ~ ■ 


T^ + A (( ^ + ^_i _ wi ,( Xl)2 + (X2)2 } 
+ «p7 « Xl)2 + (X2)2 > 2 ' 


We immediately see that the only terms from action ca involving are 

F 2“(^ L ) 2 + T ^ W ^) 2 + °( r )> 

^34 “l 2 

and therefore we can integrate them out in the limit i? —> 00 . In the leading order in R 
the remaining terms give the standard action 

5 = "2 9 202 02 f([ X l» X *] - 


for a gauge theory on Mg for general f9 //;/ . The X^ are interpreted as “covariant coordi¬ 
nates”, which can be written as 7 

Xfj, .— x^ T iO^i/Ay. 


7 We do not distinguish between upper and lower indices. 
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Hence the gauge fields A M describe the fluctuations around the vacuum. In particular, 
note that our regularization procedure clearly fixes the rank of the gauge group, unlike 
in the naive definition on Mg as discussed in [4], The generalization to the U(n) case is 
obvious. 


7.2 U( 1) Instantons on Mg 

The construction of instanton solutions for the two-dimensional noncommutative plane 
given in [4] can be easily generalized to the four-dimensional case. We shall recall and 
discuss these 4-dimensional “fluxon” solutions in some detail here, in order to understand 
the relation with the above solutions. To simplify the following formulas, we restrict our 
discussion from now on to the selfdual case 


OflV = 


'-fiispa @ pa 


and denote 


6 012 — 034j 


the generalizations to the antiselfdual and the general case are obvious. Then the action 
for {7(1) gauge theory on Mg reads 


S — 2g 2 0 2 ( 77 ) 

where 

F IM , = i([X IM ,X v ]-i0 lu ,) (78) 

is the field strength. In terms of the complex coordinates 


x± L := xi ± ix 2 , x± R := x 3 ± ix 4 , 
the commutation relations (1721) take the form 

[x +a ,x_ b \ = 2 65 ab , [; x +a ,x +b ] = [x- a ,x_ b ] = 0, (79) 

where a, b E {L, R}. The Fock-space representation 7i of (1791) has the standard basis 

\n 1 ,n 2 ), ni,n 2 G N, 

with 

X-L\ni,n 2 ) = \p2B\/n\ + l|ni + l,n 2 ), x + L\ni,n 2 ) = \J29^/iv[\ni — l,n 2 ) 
X- R \ni,n 2 ) = \flB\Jn 2 + l|ni,n 2 + 1), x +R \ni,n 2 ) = V26^n^\ni, n 2 - 1). 

Similarly, using the complex covariant coordinates X± a 

X±l = X] ± iX 2 , X± R = X%± iX\ (80) 
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and the corresponding field strength 


Xaa,/3b [^Qa: Xfjb] ^^^aP^ab 

with a, b G {L, i?} and q,/3 G {+, —}, the action (1771) can be written in the form 

2 

$ = 2Q'2 ; P'+a,—a^+a— a ~~ ^ ' X+a,+bF-a,—b)• 

y a a,b 

Then the equations of motion are given by: 

yit-Xqq, (Faa,pb)'] =0- (81) 

a,a 

Let us consider a finite dimensional subvectorspace V n of Ti of dimension n, which we can 
assume (using a unitary gauge transformation) to be spanned by a finite set of vectors 
|ni,n 2 ) G H, 

V n = {{\ik,jk)\ k = l,...,n}). (82) 

Following [4] one finds solutions to the equations of motion given by 8 

n 

X+l ■■= Sx +L S ] + ^7fcKfc, jfc)(4, jfcl (83) 

k =1 

n 

X+l ■■= Sx +R S^ + 5^7^14, jk)(ik,jk\ ■ (84) 

fc=i 

Here J^’ R G C determine the position of the fluxons, and S denotes a partial isometry 
from H to H\V n with S^S = 1, SS^ = 1 — P Vn , where 

n 

Py n := ^ 'j Nfc) jk) {’i'ki jk | 

fc=l 

is the projection operator onto the subspace V n . The held strength F fiu for this solution is 


X/J.v — PyJ^. 


In particular, the action corresponding to the instanton solution (1881841) is proportional to 
the dimension of the subspace V n 


S[4”’l = *r (P V J = ^ 


n. 


We will see in the next section that this class of solutions can be reproduced by instanton 
solutions m on Sfj L x Sjj R in a suitable scaling limit. Let us stress again that this is 
only one class of I/(l)-instanton solutions for Mg which is called “fluxons”, since they can 

8 Note that [X^l , A'$] = , X™] = [X^,X&] = [X^,X^\ = 0 . 
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be interpreted as localized flux. The localization can be seen as follows: recall [44] that 
the above projection operators can be represented on the space of commutative functions 
(using a normal-ordering prescription) as 




,X L ul ,X +L . ul .X R , k 2 X +R ,,,2 x+ L x~ L x+ R x~ R 

-(—=f (—=)''• (—==) =) e 20 20 


k l \k 2 \ y/20 V20 V20 V2 0 


Hence the above field strengths F )W = P Vn 9^ u are superpositions of Gauss-functions which 
are localized in a region in space of size \f0. 


7.3 Instantons on Mg from Sjj x S$j 

With the scaling limit of Section 17.11 the gauge theory on Sjy x S ^ provides us with a 
regularization for the gauge theory on Mg. Of course, such a regularization might affect 
the topological features of the theory, an effect we want to investigate in this section. For 
this, we will map the topologically nontrivial solutions found in Section [HI on x to 

K 

Consider again the solutions d that combine the fluxon solutions with the monopoles, 
with the fluxons at the north pole instead of the south pole because we want to study 
their structure. Their scaling limit as in d gives 


Xi 

X i+2 


[29 f diag(d^ 1 ,..., df n ) 

V N V 0 

[29 / diag (d R v ...,d R n ) 

V N V 0 


0 ) 

a L ® 1 ) ’ 

(85) 

° ^ 

a R l ® Af +m ) 

(86) 


for i = 1,2. Recalling that the rescaled Ag 2 on S'[r L x become the x-t’s on Mg in the 
scaling limit 


(X RR ± i\[ ,R ) —> x±l,r, 
we see that (1531) and become the instantons Hi on R^, 

n 

Xi+iX 2 - x[ n l = Sx +L S^ + ^2^\i k ,j k )(i k ,j k l (87) 

k=1 
n 

X 3 + iX 4 -»• = Sx +R S^ + y['Yk\ik,3k){ik,jk\- (88) 


k =1 


Here the (dj)-block acting on a basis | ik,jk) of V n C 7i = becomes the projector part 
of (1871188(1 with 


'20 ,L,R 
^20 ,L,R 

N 2 * 


T 3 L,R 

Re 7 fc , 


T L,R 

Im 7 k > 


(89) 
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and the monopole block becomes Sx+S t where S' is a partial isometry from TL to TL\V n . 
Note that we can recover any value for the 7 ’s in this scaling, solving the constraint 
didt = for d 3 Therefore the full moduli space of the fluxon solutions dHSl EH> 

on can be recovered in this way. Furthermore, the meaning of the parameters ^ L ' R is 
easy to understand in our approach: Note first that using a rotation (which acts also on 
the indices) followed by a gauge transformation, the di can be fixed to be radial at the 
north pole, d R,R ~ (0, 0,N/2). This is a fluxon localized at the north pole. Now apply a 
“translation” at the north pole, which corresponds to a suitable rotation on the sphere. 
Rotating the vector d R,R in the scaling limit amounts to a translation of the 7 ^’ R according 
to (1891) . which therefore parametrize the position of the fluxons. 

It has been noted [ 2 ] that the Sx + S t correspond to a pure (but topologically nontrivial) 
gauge, which can qualitatively be seen already in two dimensions. There, the isomorphism 
S : | k) —> |k + n) is basically (-^===) n ~ ~ e m,f and therefore the gauge field 

Ai = SdiS t has a winding number n. The topological nature of the Sx+S t is even more 
evident in our setting, as they are the limit of the monopole solutions mm on &N X $N- 
Moreover, note that their contribution to the action (EH) survives the scaling: even though 
the field strength vanishes as R —» 00 , the integral gives a finite contribution equal to the 
contribution of the fluxon part. This topological “surface term” is usually omitted in the 
literature on Mg, but becomes apparent in the regularized theory. 

So it seems that we recovered all the instantons of Section o but in fact there is an 
important detail that we haven’t discussed jet. It is the embedding of the n-dimensional 
fluxons and the (N — m)(N + m)-dimensional monopole solutions into the IV 2 -dimensional 
matrices of the ground state. Such an embedding is clearly only possible for n = m 2 . This 
means that the regularized theory has some kind of “superselection rule” for the dimension 
of the allowed instantons, a rule that did not exist in the unregularized theory 9 . 

One way to allow arbitrary instanton numbers is to allow the size A f of the matrices to 
vary. However, this is less satisfactory as it destroys the unification of topological sectors 
which is a beautiful feature of noncommutative gauge theory. On the other hand, the type 
B solutions (1531) together with the changed action © might allow the construction of 
the missing instantons. The idea is to fill up the unnecessary m 2 — n places with di = 0. 
The changed action would not suppress such solutions any more, and in fact they would 
not even contribute to the action. This amounts to adding a discrete sector to the theory 
which accommodates these type B solutions, but decouples from the rest of the model. 
Whether or not one wants to do this appears to be a matter of choice. This emphasizes 
again the importance of a careful regularization of the theory. It would be very interesting 
to see what happens in other regularizations e.g. using gauge theory on noncommutative 
tori or fuzzy CP 2 . 

9 Note that this is different in two dimensions. There, a rank n fluxon can be combined with a (N — n)- 
dimensional monopole block and all the instantons on Rg can be recovered. Furthermore, the actions for 
the fluxons and the monopoles scale differently with N. Therefore the action for the monopoles vanishes 
in the scaling limit that produces a gauge theory on Rg with rescaled coupling constant. 
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8 Fermions 

8.1 The commutative Dirac operator on S 2 x S 2 

To find a form of the commutative Dirac operator on S 2 x £ 2 which is suitable for the 
fuzzy case, one can generalize the approach of [45] for £ 2 , which is carried out in detail in 
Annendix lE.3l One can write the flat £0(6) Dirac operator D$ in 2 different forms, using 
spherical coordinates of the spheres and also using the usual flat Euclidean coordinates. 
Then one can relate Dq with the curved four-dimensional Dirac operator D± on £ 2 x £ 2 
in the same spherical coordinates. This leads to an explicit expression for O 4 involving 
only the angular momentum generators, which is easy to generalize to the fuzzy case. The 
result is rather obvious and easy to guess: 

d 4 = J J ) + i ( M = r% + 2C 0 , (90) 

which is clearly a £0(3) x £0(3)-covariant Hermitian first-oder differential operator. Here 
T^ generate the £0(6) Clifford algebra (1401) . Co is defined in (FHll . and we put R = 1 for 
simplicity here. However this Dirac operator is reducible, acting on 8 -dimensional spinors 
Tg corresponding to the £0(6) Clifford algebra. Hence Tg should be a combination of 
two independent 4-component Dirac spinors on the 4-dimensional space £ 2 x £ 2 . To see 
this, we will construct explicit projectors projecting onto these 4-dimensional spinors, and 
identify the appropriate 4-dimensional chirality operators. This will provide us with the 
desired physical Dirac or Weyl fermions. 

8.1.1 Chirality and projections for the spinors 

There are 3 obvious operators which anti-commute with D 4 . One is the usual 6 -dinrensional 
chirality operator 

r := *TfT£iirMrf = ( “ 1 J ) , ( 9 i) 

which satisfies 

{d 4 ,t} = o, r t = r, r 2 = i. (92) 

The 8 -component spinors 'hg split accordingly into two 4-component spinors 4*g = ( - 1 J, 

which transform as (4) resp. (4) under so( 6 ) = sit(4); recall the related discussion in Sec¬ 
tion QJ The other operators of interest are 

XL = T lL XiL and \R = T lR XiR- 

They preserve £0(3) x £0(3) C £0(6), and satisfy 

{D4,Xl,r} = 0 = {xl,Xr} 

as well as 

Xl,r = 1- 
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We will also use 


( 93 ) 


x = -j= r % = -j= (xl + xr) 

which satisfies similar relations. This means that 

P± = ^(1 ±iXLXR) (94) 

with 

P± = P±, P + + P- = 1 and P+P- = 0 (95) 

are Hermitian projectors commuting with the Dirac operator on S 2 x S 2 as well as with 

r, 

P± = P± and [P±, D 4 ] = [P±,T] = 0. (96) 

Therefore they project onto subspaces which are preserved by D 4 and T, and are invariant 
under SO( 3) x SO (3). Hence the spinor Lagrangian can be written as 

= T 1 f D 4 ^ + + D 4 T_ 

involving two Dirac spinors 'h-t = P±'$>s. In order to get one 4-component Dirac spinor, 
we can e.g. impose the constraint 

P + T 8 = \H 8 , (97) 

or equivalently give one of the two components a large mass, by adding a term 

8 (98) 

to the action with M_ —> 00 . The physical chirality operator is now identified using m 
and m as T acting on f + . It can be used to define 2-component Weyl spinors on S 2 x S 2 . 
To make the above more explicit, consider the north-pole of the spheres, i.e. 



and xr = 



In the basis m for the Clifford algebra we then get explicitly 


P ± = 2 (1± * 


" 7 l 7 r 


j ) = ^(1 ±< r 3 ^ 3 )- 

ILIR J 1 


This means that 

P+ = diag(l, 0,0,1,0,1,1,0) 

projects onto a 4-dimensional subspace exactly as expected. 
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8.2 Gauged fuzzy Dirac and chirality operators 

To find fuzzy analogues of (Em and ESI coupled to the gauge fields, we recall the connec¬ 
tion between the gauge theory on S 2 x S 2 and the 50(6) Gamma matrices established in 
Section In the spirit of that section a natural fuzzy spinor action would involve 

&CV, ( 99 ) 


where T is now a 8A f x jV-matrix (with Grassman entries). Of course, (1991) does not have 
the appropriate commutative limit, but we can split C into a fuzzy Dirac operator D and 
the operator y defined by 

= ( 100 ) 

which generalizes ESI); we used here the definition E2J),(ES1) of Co- This operator satisfies 



and reduces to ESI) in the commutative limit. Note also that y commutes with gauge 
transformations, since the coordinates A^ are acting from the right in mm- Setting 

J> = [A m ,T], 


we get for the fuzzy Dirac operator 


D = C-^=x = T^J, + A,) 


+ 2 Co = + 2 Co. 


( 101 ) 


Here 10 


Jn + 


( 102 ) 


is a covariant derivative operator, i.e. = Ty^Uip which is easily verified using ESI) . 

This D clearly has the correct classical limit m for vanishing A. and the gauge fields are 
coupled correctly. In particular, this definition of D applies also to the topologically non¬ 
trivial solutions of Section El without any modifications. Moreover, the chirality operator 
T as defined in ED anti-commutes with D also in the fuzzy case, 


{D,T} = 0. (103) 

In particular there is no need to consider e.g. fuzzy Ginsparg-Wilson operators as in the 
2-dimensional case [46-48]. However, the anticommutator of D and y no longer vanishes. 
We find 

0, X} = ~^§ (2(A /i + H M )J /i -2H /i A /i + {r^,C' 0 }P M + 2) =0(1), (104) 

10 We set R = 1 in this section for simplicity. 
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since x^J^, = O(i) and = O(^) using (128f) . Furthermore, using some identities 

given at the beginning of Section 0 we obtain for D 2 i/j: 

D 2 il> = {^F^ + V fl V lx + {T^C 0 }V fl + 2)i; 

= : Ffv, + □ + 2)V>, (105) 

defining the covariant 4-dimensional Laplacian □ acting on the spinors. This corresponds 
to the usual expression for D 2 on curved spaces, and the constant 2 is due to the curvature 
scalar. Since D 2 and Yi ,w F fW are both Hermitian and commute with T and P± as defined 
in dmzD in the large N limit, it follows that □ satisfies these properties as well. Note that 
HMD can also be written as 

(-D - C 0 ) 2 = V u 'F, lv + V t ,V, L + 1 (106) 

which might suggest to interpret T> II T> II as covariant Laplacian; however this is not correct 
since T>^D^ does not commute with the projections P± HMD even in the commutative 
limit. The reason for this is our formulation using spinors based on the 50(6) Clifford 
algebra rather than 50(4) spinors and comoving frames. The corresponding projections 
to physical Dirac- or Weyl-spinors in the fuzzy case will be discussed next. 

8.2.1 Projections for the fuzzy spinors 

For the fuzzy case, we can again consider the following projection operators 

x.l* = ^(r lL TA lL + C 0 L T), 
xr* = l(r ii? vi/A ii? + o ( fT) 

which satisfy 

Xl,r = 1 , {xl,Xr} = 0. 

This implies ( XlXR ) 2 = ~ L and we can write down the following projection operators 

^± = ^(i±*xafi) (107) 

which have the classical limit (IMll and the properties However, the projector no 

longer commutes with the fuzzy Dirac operator Gnu: 

0 , XLXR } = {D, Xl}xr - Xl{D, Xr} 

= — jy + AiL)JiL ~ 241 h\l + 2Cq j T lL T>iL + 1 )XR 

—Xl(2(\r + AiR)JiR — 2Aiu\iR + 2C§ T lR T>iR + 1 )^ , 

which only vanishes for N —► oo and tangential A^ To reduce the degrees of freedom 

to one Dirac 4-spinor, we should therefore add a mass term 

(108) 
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which for M_ —> oo suppresses one of the spinors, rather than impose an exact constraint 
as in (1371). This is gauge invariant since P± commutes with gauge transformations, 

P±ip^ UP±ip. 

The complete action for a Dirac fermion on fuzzy x Sfj is therefore given by 

SDirac = J (109) 

with M_ — * oo. The physical chirality operator is given by rm, which allows to consider 
Weyl spinors as well. 

9 Conclusion and outlook 

We have constructed U(n) gauge theory on fuzzy S jy x as a multi-matrix model. 
The model is completely finite, and can be considered as a regularization either of Yang- 
Mills on the commutative S 2 x S 2 , or on the noncommutative Mg in a suitable scaling 
limit. The quantization is defined by a finite “path” integral over the matrix degrees of 
freedom, which is convergent due to the constraint term. A gauge-fixed action with BRST 
symmetry is also provided. We then discussed some topologically non-trivial solutions in 
the 17(1) case, which reduce to the known “fluxon” solutions on Mg in the appropriate 
scaling limit, reproducing the full moduli space. On x Sfj they arise as localized 
flux tubes together with a monopole background field. This provides a very clean non- 
perturbative definition of noncommutative gauge theory with fixed rank of the gauge 
group U(n), and a simple description of instantons as solutions of the equation of motion 
in one single configuration space. Furthermore, we have shown how charged fermions in 
the fundamental representation can be coupled to the gauge field, by defining a suitable 
Dirac operator D. This is easily extended to Weyl fermions using a chirality operator 
which exactly anticommutes with D. All this supports the programme to formulate and 
study physically interesting models on noncommutative spaces. 

There are many interesting conclusions and applications to be explored. One crucial 
feature is the fact that the model is completely regularized, i.e. the quantization is well- 
defined without any divergences for finite N. This should allow to study suitable scaling 
limits in N in a well-defined framework, and the emergence of an interesting low-energy 
limit which could be either commutative or noncommutative. Such a matrix regularization 
is very interesting in view of the UV/IR mixing, which indicates a close relationship be¬ 
tween NC field theory and matrix models. For example, one might try to extend the results 
in [49] in this context. We also explored some alternative formulations using “collective 
matrices” based on SO(6). Such formulations are possible only in the noncommutative 
case, and lead to the hope that new non-perturbative techniques in the spirit of random 
matrix theory may be developed along these lines. 

Another important aspect is the coupling to fermions, which could be extended to 
scalars and allows to study spontaneous symmetry breaking and the possible generation 
of other gauge groups in the low-energy limit. Finally, a detailed comparison with other 
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finite models of NC gauge theory in 4 dimensions such as [22,24,25] would be very desirable, 
to see which features are generic and which are model-dependent. 
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A The standard representation of the fuzzy sphere 

The irreducible A-dimensional representation of the su( 2) algebra Aj Q is given by 

(A 3 )« = S kl N + 1 2 ~ 2k , (110) 

(A +)ki = bk+i,i\J (A - k)k, (111) 

where k, l = 1,..., N and A± = Ai ± 1 X 2 - 

B Alternative formulation using 4J\T x 4J\T matrices 

Let us rewrite the action (SZl) in terms of the AN x AN matrices Bl, B r Noting that 
CiC R + C R C t = ( (112) 

we can rewrite Sq (H51) as 

/ N 2 \ 2 

S 6 = 2Tr (b 2 l -B 2 r - — \ + 2Tr ([B L , B R } 2 ) , (113) 

where the trace is now over AN x AN matrices. Similarly 

S break = - 2 TV (s£ - (-B% - ^ (114) 

and combined we recover m as 

/ N 2 N 2 \ 

5 = 5 6 - 25 break = 2TV l(B 2 - — f + ( -B 2 R - —f + [B L , B R } 2 \ . (115) 

This looks like a 2-matrix model, however the degrees of freedom Bl,B r are still very 
much constrained and span only a small subspace of the AN x AN matrices. We would 
like to find an intrinsic characterization without using the 7 ^ explicitly. One possibility is 
to choose the 7 ^ to be completely anti-symmetric matrices, see Appendix^! However this 
does not extend to B, since the B^ should be Hermitian and not necessarily symmetric, and 
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moreover the 7 ^ are not Hermitian (the conjugate being the intertwiner ( 6 ) C (4) < 8 ) (4)). 
Another possibility is provided by the following representation of the 7 -matrices: 

7 l = o'*® hx2, 1 r = hx 2 < 8 ) icr 1 . (116) 

They satisfy the relations pi) - (32), but are not antisymmetric. Now note that 

7 l R = iPl l L P (117) 

where 

/ 1 0 0 0 \ 

0 0 i 0 

0 10 0 

\ 0 0 0 1 / 

permutes the two tensor factors and satisfies 

P 2 = 1. (119) 

Therefore we can characterize the degrees of freedom in terms of 2 Hermitian 2J\[ x 2J\[ 
matrices 

X L = BicJi + 1 X R = B\ t a % + 1 (120) 

which are arbitrary up to the constraint that X^ L n = \- Then 

Bl = Xl®1 2 x 2 , Br = iP(XR®\<2x2)P\ ( 121 ) 

they could be extracted from a single complex matrix B = (Xl + {Xr) ® 12 x 2 - Further¬ 
more, matrices of the form X (g> 12x2 are characterized through their spectrum, which is 
doubly degenerate; indeed any such Hermitian matrix can be cast into the above form 
using suitable unitary SU(4Xf) transformations. Similarly, P can also be characterized 
intrinsically: any matrix P written as 


= 5(1 + 


a 


(118) 


P = Pq < 8 > 12x2 + Pi ® <y l (122) 

which satisfies the constraints 

Po = P 2 = 1 (123) 

is given by ceiej up to an irrelevant unitary transformation U <g> 1. We could therefore 
write down the action lira in terms of three matrices Br, —IPBrP and P, all of which 
are characterized by their spectrum and constraints of the form (..)o = The hope is 
that such a reformulation may allow to apply some of the powerful methods from random 
matrix theory, in the spirit of [7]. However we will leave this for future investigations. 
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C Stability analysis of the SO (6) - invariant action (14511 


Consider the action (1351) . We will split off the radial degrees of freedom for large N by 
setting R = 1 and 11 

Bn = A n + AiL = A iL + An + x iL $ L 

requiring that A hAil = 0, and similarly for Bm, The stability of our geometry will depend 
on the behavior of d> L and d> R . We calculate that 

.V 2 - 1 1 

Bfj.Bfj, -—— = + d>^) + dqd>£ + d>/jd>^ + Af+Afj, — [A m , A + 0(—), 

where we used that A i a Ai a = 0 and therefore both Ai a Xi a = O( jj) and Ai a [\i a , •] = O(jj) 
for a = L, R. Setting 


d>L + $R = $ 1 , 

~ ®Ft = d > 2 


we get 

IV 2 — 1 1 

B^B^ - —— = IV dq + dqdq + dqdq + A/jAn — [Xu, A + O(-). (124) 

In the limit IV —> oo we can integrate out dq, as it acquires an infinite mass. Alternatively 
we can rescale dq by setting </>i = -^dq. Then, all the terms involving </>i but the first one 
in will be of order and we can equally integrate out cf) i- 

The terms from 

FijF/i - [Bn,B iR ] 2 

involving the remaining dq will be (in the limit IV —* oo) 

-d> 2 d >2 — ^i(d> 2 )J /i (d> 2 ) - FiLXiL®2 + FiRXiR^-2 
with the tangential derivatives Jj a = — RijkXjadka■ Calculating that 

>//rdq J^dq = — c^dqc^dq — £iL$jLdqx.,x<9jLdq — XiRdiR&2XjRdjR<&2 
and using partial integration under the integral this gives 

-dqdq - dqc^< 9 /t dq — x t j J d t j j d> 2 xjj J d-j^ d > 2 — XiRdiR&2XjRdjR&2 — FiLXn<&2 + FiRXiR<&2 

Expanding both dq and F in left and right spherical harmonics as 

$2 =Y1 C klmn Y k m Y ln and F la X la = ^ fklmn Y km Y ln 
klmn klmn 

n The fact that this leads to non-hermitian fields for finite N is not essential here. 
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we get for fixed klmn, setting c = c^ mn , /“ = f%i mn and p = \ + l (l + 1 ) + k(k + 1 ) the 
following expression 

pc 2 - cf L + cf R = p(c - ±f L + Z-f R f f R ?. 

Zp Zp 4 p 

Integrating out the c’s and putting everything back this leaves us with the additional term 


-( FiL'XiL - F iR XiR)—j -—— 

4(^ - df.d,. 


(F iL XiL - F iR XiR ) 


in the action m- 


D Representation of the 50 ( 6 )- intertwiners and Clifford 
algebra 


We will use the Pauli matrices 


a 1 = 


0 1 
1 0 


u 2 = 


0 -i 
i 0 


u 3 = 


1 0 
0 -1 


which satisfy 


oV = S ij + ie ijk a k . 

With these we define the 4-dimensional antisymmetric matrices 

= CF 1 ® (7 2 , 7£ = fJ 2 (g)l, 7^ = u 3 ® cr 2 , 

7 ^, = i a 2 <g) a 1 , ^ = i 1 <g> a 2 , 7 ^ = i a 2 <g) a 3 . 


(125) 


(126) 


They are the intertwiners between SU (&) ® SU(4) and SO(6 ) and fulfill the following 
relations: 


hi)* = 7I1 
(Tk) f = - 


~'R 


and 


iH = 

7 Wr = ?t£, 

[7i; Tr] = 0 . 

We can now define the 8 -dimensional representation of the S'0(6)-Clifford algebra as 

0 Y 

7^t 0 

with the desired anticommutation relations 

7^7 v t + 7^7^! 0 


r /'< _ 


( 127 ) 


(r^, r^} = 


0 




= 2<W7 
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The chirality operator in this basis is 

r — ,'r 1 r 2 r 3 r 1 r 2 r 3 — ( ^ ^ ^ 

1 — 11 L l L l L L R L R l R— l Q J !• 

The 8-dimensional S'0(6)-rotations are generated by 

y nv = _h vt i v v} = _!_( 7 / V t - 0 \ 

8 41 1 ’ 1 1 4 l 0 7 / 4 7 " _ ydy* }■ 

E The Dirac operator in spherical coordinates 

For a general Riemannian manifold with metric 

g = g l u,dx' l dx 1 ' 

the Christoffel symbols are given by 

= \g cX {dng\ 1/ + - d \9nv)- (128) 

We can change to a non-coordinate basis (labeled by latin indices in contrast to the greek 
indices for the coordinates) by introducing the vielbeins e« with 

g„v = e“e% b , g^ = e^ b 5 ab . 

With these, the Dirac operator is given by 

D = -i7 a e£(«9 M + ^ a6 [ 7 a ,7 6 ]), 
where the 7 “ form a flat Clifford algebra, i. e. 

{ 7 “,7 6 } = 2 5 ab , 7 at = 7 a 
and the spin connection oj fulfills 

^“-^ + <6^ = 0. (129) 

E.l The Dirac operator on M 6 in spherical coordinates 

We will now write down the flat SO (6) Dirac operator Dq by splitting M 6 into x 
and introducing spherical coordinates on both the left and right hand side. The flat metric 
becomes 

= r\ dd L <g> d0 L + r 2 sin 2 d L d(j>L ® #l + dr L ® dr L 

+ r R dOR <g> d0R + rf { sin 2 Or d(j)R <g> d<f>R + drR ® drR. 


g e 
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(130) 

(131) 



Looking at the formula for the Christoffel symbols mm we see that all the symbols with 
both right and left indices vanish. For the symbols with only right or only left indices we 
get 


•n 

-e- 

II 

— sin 9 cos 9, 

(132) 

p '/ 1 _ 

1 e<j> ~ 

cos 9 0 

sin 9 ~ * e ’ 

(133) 

r r — 

1 ee ~ 

-r, 

(134) 

r r — 

1 H 

—r sin 2 9, 

(135) 

II 

05 

05 S- 

[-H 

~j i h-> 

II 

Tl 

■i 

(136) 

p'A _ 

1 r<j> ~ 

^ _ p</> 
r </>»•> 

(137) 


where we have dropped the left or right subscript for simplicity. All other symbols vanish. 


We want to go to a non-coordinate basis by introducing the vielbeins 

= r L\ e l L L = r L sin 9l\ eh = ^ ( 138 ) 

4n = r l; = r R sm0 R - = 1. (139) 

Calculating the spinor connection by cm we again see that all the terms with both left 

and right indices vanish. The terms with only left or only right indices are 

^0 2 = -cos0 =-w / l5 (140) 

w 0 2 3 = sin 6 =-w/ 2 , (141) 

^0*3 = 1 = — (142) 


where we again dropped the left or right subscripts. Putting all this together we see that 
D$ splits up into a left part D 3 l and a right part D 3R as 

D§ = D 3 l + D 3R (143) 


with 


D 3 l 

d 3R 


■u 1 1 la , cos ^L -r 2 1 q 

— — \do L + ■ /) ) ~ -:— 7T^0l 

pl sm.0L 'f'h smtT 


- + P) 


-ir 


i i 


(^h + 


cos 0 ft , 


R r R y "° R ' sin 0 ft 


-ir 


1 




r/? sm 0 R VR r R 


where the T have to form a SO(Q ) Clifford algebra. 


(144) 

(145) 


E.2 The Dirac operator on S' 2 x S 2 

We now want to calculate the curved Dirac operator on S 2 x S 2 in the spherical 
coordinates of the spheres (they are the same spherical coordinates we used before, now 
restricted to the spheres). The metric on S 2 x S 2 with radii rp and r R is 

g A = r- 2 d&L ® d0L + r 2 L sin 2 9 l d(j)L < 8 > d^L 

+r R d0 R <g> d0 R + r R sin 2 0 R d<j) R <g> d(f) R . 
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The metric is the same as (IT3U1) restricted to the spheres, so the Christoffel symbols are 
the same as (Usui and dim Again introducing the vielbeins 

e\ L L =r L ] = r L sin 0 L - (146) 

ej" = r L ; e = r R smO R , (147) 

we see that also the spin connection is the same as CHI, and therefore we can again split 
D 4 into a right part D 2 R and a left part D 2 L as D 4 = D 2 L + S> 2 R with 


D 2 L 

D 211 


- ir L — ( d ° L + 

-iT l R —{de R + 
r R 


cos Ql 
sin 01 
cos 6 r ^ 
sin 0 R 


- ?T 


1 


sin 

^<pR, 


— ?T 2 
u R 


r R sin 9 r 


where the T form a flat 50(4) Clifford algebra. 


(148) 

(149) 


E.3 50(3) x SO (3)-covariant form of the Dirac operator on 5 2 x 5 2 

The flat 50(6) Dirac operator Dq was split into a left part D 3R and a right part D 3R using 
spherical coordinates in dm . Of course, Dq can also be written in the usual Euclidean 
coordinates as 

Dq = -zT^, 

where again we can split it into a left and a right part as 

Dq = D 3 l + D 3R 


with 


D 3 l = —iT l L di, D 3R = — iO R di, (150) 

{D 3 l,D 3R } = 0. (151) 

We have left open which representation of the 50(6) Clifford algebra we want to use for 
the T in (II4411 4511 . but T in csii is really the representation given by m- We will now 
relate the two expressions for the Clifford algebra and the Dirac operator by first defining 


JiL — S.jjkXj i/^kL and Jt R — i^ijkXj Rdk R 


and noting that 


We calculate that 


2 

r L d iL 



r l R d iR 


rwy 1 

rL J V r R J 


( T L X 3l \ f XildiL _ }_ f Yl 0 
\ r L J \ r L rt, \ 0 7l 



( \ ( Xi R d iR J_( Yr 0 

\ r R J \ r r r R \ 0 7 ^ 



(152) 


(153) 

(154) 


34 
















and therefore 


D3L 


Dsr 


—i 


T L X 3L 

tl 



1 

tl 



7 l L 

0 


7 5* 



(155) 

(156) 


Comparing this with 114411451) we see that 



and 





(157) 


as the Jl and Jr have no radial components. From (I155I156[) we can also deduce that 


I )]=0=[ r K,( _°1 (158) 

and 

r H ,(; ;) } =°= { n,(_° 1 <»«» 

The curved Dirac operator D 4 on S 2 x 5 2 expressed in the spherical coordinates of the 
spheres also splits up as D 4 = D 2 L + D 2 R with right part D 2 R and left part D 2 L given in 
cm , cm - Comparing this with 114411451) , we see that the dependence on the tangential 
coordinates is the same in both expressions. With (115811591) we see that the matrices 

^ ^ t\t 1 l and ^ ^ q ^ for i.j = 1,2 form a 50(4) Clifford algebra and 

can therefore be used as the T. Note that this representation is still reducible, a problem 
we deal with in Section f 8 . 1.11 Now we can get a simple relation between the O 3 restricted 
on the spheres and the D 2 



“(l 0 ) (* r ^3LU.-^) = D 2L , 

-i ^ q ^ (irRD3R\res. ~ ~) = D 2 R. 

Inserting (11551156(1 and using HSU together with (fl52l) we find that 

D 2 l = J J)), (160) 

D 2 r = ^(r Ufl + i (_° 1 J))- ( 161 ) 

Setting ri = tr = 1 for simplicity, the Dirac operator O 4 on 5 2 x 5 2 takes the form Pip . 
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